净 多 样 化 速率 和 进化 时 间 对 虎 耳 草 目 科 间 物 种 多 样 性 差异 
的 影响 
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摘 要 : 不 同 生物 类 群 包 含 的 物种 数目 常 存在 巨大 差异 ， 这 是 生态 学 和 生物 学 研究 中 普遍 
观察 到 的 现象 。 然 而 ， 这 一 现象 产生 的 原因 仍然 是 未 解 之 恋 。 从 宏观 进化 的 角度 ， 进 化 时 


间 假 说 和 多 样 化 速率 假说 是 两 个 比较 流行 的 假说 。 进 化 时 间 假 说 认为 类 群 的 演化 时 间 越 长 ， 
积累 的 物种 丰富 度 越 高 ， 而 多 样 化 速率 假说 认为 类 群 的 净 多 样 化 速率 越 快 ， 则 其 物种 丰富 
度 越 高 。 为 验证 这 两 个 假说 ， 该 文 以 一 棵 包含 1 539 个 物种 化 石 定年 的 虎 耳 草 目 系 统 发 育 
树 为 基础 ， 通 过 宏观 进化 分 析 获 取 了 虎 耳 草 目 内 15 个 科 的 物种 形成 和 灭绝 速率 ， 并 计算 了 
每 个 科 的 平均 多 样 化 速率 。 结 果 表 明 : (1) 虎 耳 草 目 的 物种 多 样 化 速率 有 着 增加 的 趋势 ， 并 
且 多 样 化 速率 的 增加 主要 出 现在 温带 和 高 山 类 群 ， 如 茶 若 子 科 、 景 天 科 和 芍药 科 等 。(2) 采 
用 系统 发 育 广 义 最 小 二 乘 模型 和 线性 回归 模型 结果 表明 ， 虎 耳 草 目 15 个 科 的 物种 丰富 度 与 
科 的 分 化 时 间 和 科 内 物种 的 最 近 共 同 祖先 年 龄 都 没有 显著 相关 关系 ， 而 与 净 多 样 化 速率 显 
著 正 相关 (R =0.38，P<0.01)。 该 研究 支持 了 多 样 化 速率 假说 ， 认 为 不 同 科 的 净 多 样 化 速率 
的 差异 是 导致 虎 耳 草 目 科 间 物种 数目 差异 的 主要 原因 之 一 。 全 球 气候 变 冷 可 能 为 虎 耳 草 目 
中 草本 、 落 叶 乔 木 和 灌木 等 能 够 适应 寒冷 环境 的 类 群 提供 了 分 布 范围 扩张 和 物种 快速 多 样 
化 的 机 会 。 该 研究 表明 在 温带 和 高 山 扩张 类 群 中 ， 物 种 净 多 样 化 速率 可 能 是 导致 不 同类 群 
物种 数目 差异 的 主要 原因 。 

关键 词 : 虎 耳 草 目 ， 物 种 多 样 性 ， 系 统 发 育 ， 多 样 化 速率 ， 宏 观 演化 
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greatly, which is a phenomenon commonly observed in ecological and biological studies. However, 
the underlying mechanisms of this pattern are still in hot debate. From the perspective of 
macroevolution, time-for speciation hypothesis and diversification rate hypothesis are two popular 
hypotheses. Time-for speciation hypothesis suggesting that the time-for-speciation is the primary 
driver of diversity variation among groups by species accumulation along time, while 
diversification rate hypothesis suggesting that the groups with high species number resulted from 
their high net diversification rate. To test these two hypotheses, this research studied species 
variation among 15 families of Saxifragales based on a fossil-calibrated phylogenetic tree 
containing 1 539 species. We obtained the rates of speciation and extinction of 15 families and 
calculated the average rate of diversification of each family through macro-evolutionary analysis. 
The results are as follows: (1) The species diversification rate of the order Saxifragales shows an 
increasing trend, and the increase of diversification rate mainly occurred in temperate and alpine 
groups, such as Grossulariaceae, Crassulaceae and Paeoniaceae. (2) The results of phylogenetic 
generalized least squares and linear regression model show that the species richness of the 15 
families is not significantly correlated with the age of family differentiation or the age of the most 
recent common ancestor, but it is positively correlated with the rate of net diversification (R°=0.38, 
P«0.01). These results support the diversification rate hypothesis, suggesting that the 
diversification rate is a primary evolutionary driver of the variation of species richness among 
families within Saxifragales. Global cooling may provide opportunities for the expansion of 
distribution and rapid species diversification of herbaceous, deciduous trees and shrubs that can 
adapt to cold conditions. The study showed that in temperate and alpine spreading groups, the net 
diversification rate may be the main reason for the difference in species numbers. 


Key words: Saxifragales, species diversity, phylogeny, diversification rate, macro evolution 


生物 类 群 之 间 的 物种 总 数 的 差异 是 自然 界 中 普遍 存在 的 现象 ， 但 是 这 一 现象 的 形成 机 
制 仍然 是 未 解 之 谜 (Li & Wiens, 2019)。 从 生物 类 群 的 宏观 进化 角度 ， 进 化 时 间 假 说 (time-for 
speciation hypothesis) 和 多 样 化 速率 假说 (diversification rate hypothesis) 是 解释 类 群 间 的 物种 
多 样 性 差异 的 两 个 最 为 核心 的 假说 (Scholl & Wiens, 2016; Li & Wiens, 2019)。 进 化 时 间 假 说 
认为 ， 进 化 时 间 越 长 ， 物 种 能 积累 的 物种 多 样 性 越 高 ， 因 此 进化 的 时 间 差 异 是 导致 类 群 间 
物种 多 样 性 差异 的 主要 原因 。 换 名 话说， 一 个 起 源 古 老 的 类 群 ， 比 新 近 起 源 的 类 群 拥有 更 
长 的 进化 时 间 ， 因 此 可 以 积累 更 多 的 物种 。 多 样 化 速率 假说 认为 ， 在 生物 演化 历史 中 ， 不 
同类 群 种 化 速率 和 灭绝 速率 的 动态 变化 导致 的 净 多 样 化 速率 的 差异 是 引起 物种 多 样 性 变化 
的 主要 原因 (Scholl & Wiens, 2016)。 也 就 是 说 ， 物 种 形成 速率 高 而 灭绝 速率 低 的 类 群 具 有 更 
高 的 净 多 样 化 速率 ， 因 此 具有 更 多 的 物种 。 
目前 ， 进 化 时 间 假 说 和 多 样 化 速率 假说 在 物种 多 样 性 形成 中 的 相对 作用 还 存在 很 多 争 
议 。 进 化 时 间 假 说 成 立 的 前 提 是 演化 过 程 中 净 多 样 化 速率 在 不 同类 群 之 间 没 有 差异 或 者 差 
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异 极 小 。 然 而 越 来 越 多 的 研究 发 现 ， 净 多 样 化 速率 在 不 同类 群 之 间 有 着 差异 。 特 别 是 经 历 
过 大 型 灭绝 事件 的 类 群 ， 物 种 数目 会 又 然 下 降 ， 从 而 导致 物种 多 样 性 的 突然 丧失 ， 物 种 数 
目 减少 。 例 如 ，Rabosky 等 (2012) 对 多 细胞 真 核 生 物 1 397 个 主要 分 支 的 研究 发 现 ， 分 支 的 
分 化 时 间 ， 即 干 龄 (stem age)， 与 分 文 所 包含 的 物种 数 之 间 没 有 显著 关系 ， 因 而 否定 了 进化 
时 间 假 说 。Pyron & Wiens (2013) 发 现 ， 许 多 古老 两 栖 动 物 分 支 具有 较 低 的 物种 多 样 性 ， 但 
是 一 些 年 轻 的 分 支 具有 较 多 的 物种 。 但 是 该 研究 并 没有 直接 验证 物种 多 样 性 与 多 样 化 速率 
之 间 的 关系 。 而 McPeek & Brown (2007) 对 不 同 动物 类 群 物种 多 样 性 的 研究 却 发 现 分 支 年 龄 
是 物种 多 样 性 的 决定 因素 。Marin & Hedges (2016) 对 不 同类 群 的 冠 龄 (crown age) 与 物种 丰富 
度 的 关系 的 研究 也 表明 两 栖 动 物 、 鸟 类 和 哺乳 动物 的 物种 丰富 度 与 其 最 近 共 同 祖先 的 年 龄 
(age) 显 著 正 相关 。 为 进一步 验证 进化 时 间 和 多 样 化 速率 假说 ，Scholl & Wiens (2016) 使 用 生 
命 之 树 (tree of life) 对 不 同 生 物 类 群 不 同 分 类 学 等 级 的 物种 多 样 性 的 差异 进行 了 研究 。 研 究 
中 发 现 ， 在 不 同 的 分 类 等 级 ， 例 如 ， 门 、 目 和 科 等 ， 物 种 多 样 性 主要 是 由 多 样 化 速率 决定 
的 ， 而 与 进化 时 间 没 有 显著 相关 关系 ， 甚 至 出 现 负 相关 关系 。 古 老 的 类 群 物种 多 样 化 速率 
低 ， 而 年 轻 类 群 的 多 样 化 速率 高 可 能 是 导致 该 现象 的 主要 原因 。 并 且 ， 进 化 时 间 与 物种 多 
样 化 速率 的 关系 在 不 同 分 类 等 级 和 不 同 的 生物 类 群 中 的 差异 可 能 导致 不 显著 的 物种 多 样 化 
速率 与 物种 多 样 性 的 关系 。 综 合 来 说 ， 进 化 时 间 假 说 和 多 样 化 速率 假说 在 物种 多 样 性 形成 
中 的 作用 仍然 存在 争议 。 

虎 耳 草 目 (Saxifragales) 包 含 了 15 个 科 的 共 约 3 000 个 种 ， 包 括 了 和 乔木， 灌木 ， 多 年 生 
> 或 一 年 生 草 本 植物 ， 多 肉 植物 和 水 生 植物 等 多 种 生活 型 ， 广 泛 分 布 于 全 球 不 同 的 生态 系统 
Se (Soltis et al, 2013)。 虎 耳 草 目 内 部 科 的 分 化 时 间 从 早 和 白垩 纪 晚期 至 第 三 纪 早 期 ， 且 科 内 
物种 数目 存在 巨大 差异 ， 最 大 的 景 天 科 包 含 了 大 约 1 600 个 种 ， 而 锁 阳 科 (Cynomoriaceae)、 
PY .t» ACF} (Tetracarpaeaceae) All ba Ht RE} (A phanopetalaceae)A] KB AH 10 个 物种 。 同 时 ， 虎 
耳 草 目 的 分 子 系统 学 研究 基础 扎实 ， 基 于 301 个 核 基 因 建 立 的 系统 发 育 骨 架 ， 确 定 了 各 科 
之 间 的 系统 发 育 关 系 和 分 化 时 间 (Eolk et al., 2019). Folk 等 (2019) 还 结合 小 片段 构建 了 种 级 
水 平 的 系统 发 育 树 ， 物 种 覆盖 率 达 到 当时 公认 物种 70% 以 上 ， 为 准确 估计 物种 的 多 样 化 速 
率 黄 定 了 基础 。 因 此 ， 虎 耳 草 目 是 研究 科 间 物种 丰富 度 差 异 、 验 证 进化 时 间 和 多 样 化 速率 
假说 的 理想 类 群 。 本 研究 采用 系统 发 育 广义 最 小 二 乘 模 型 (phylogenetic generalized least 
squares，PGLS) 分 析 了 虎 耳 草 目 科 的 物种 多 样 性 与 科 的 冠 龄 、 干 龄 和 多 样 化 速率 的 关系 ， 
发 现 了 物种 多 样 性 与 多 样 化 速率 有 显著 的 正 相 关 关 系 ， 与 进化 时 间 没 有 显著 关系 ， 且 进化 
时 间 与 多 样 化 速率 的 综合 效应 对 这 种 物种 丰富 度 差异 有 更 强 的 解释 性 。 该 结果 发 现 了 多 样 
化 速率 在 虎 耳 草 目 物种 多 样 性 形成 中 的 重要 作用 ， 支 持 了 物种 多 样 性 格局 的 多 样 化 速率 假 
说 ， 同 时 也 支持 了 进化 时 间 和 多 样 化 速率 的 综合 效应 。 

1 材料 与 方法 
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1.1 虎 耳 草 目 系统 发 育 树 
本 研究 采用 的 系统 发 育 树 主要 基于 Folk 等 (2019) 年 发 表 的 虎 耳 草 目的 系统 发 育 树 和 徐 
莹 等 (2021) 年 构建 的 虎 耳 草 属 的 系统 发 育 树 。Folk 等 (2019) 首 先 基于 来 自 627 个 物种 的 301 
个 核 基因 序列 ， 构 建 了 虎 耳 草 目 主要 分 支 的 系统 发 育 树 并 通过 化 石和 分 子 钟 定年 确定 了 主 
要 分 支 的 分 化 时 间 。 之 后 ， 基 于 GenBank 中 下 载 虎 耳 草 目的 核 基因 ITS 片段 和 叶绿体 基因 
组 的 matK 等 24 个 基因 片段 ， 建 立 了 包含 15 个 科 1 455 个 物种 的 系统 发 育 树 。 徐 莹 等 (2021) 
构建 的 系统 发 育 树 包含 了 353 个 虎 耳 草 属 的 物种 ， 涵 盖 了 Folk 等 (2019) 构 建 的 虎 耳 草 目 系 
统 发 育 树 中 的 所 有 虎 耳 草 属 物种 ， 虎 耳 草 属 的 分 化 时 间 也 与 Folk 等 (2019) 一 致 。 因 此 ， 本 
研究 将 两 棵 系统 发 育 树 进行 整合 ， 在 R 语言 中 将 徐 莹 等 (2021) 构 建 的 虎 耳 草 属 系统 发 育 树 
替换 了 Folk 等 (2019) 虎 耳 草 目 系 统 发 育 树 中 虎 耳 草 属 分 支 (Bininda & Olaf, 2004; Gordon, 
1986)。 最 终 整 合 得 到 了 一 棵 包含 1539 个 种 的 虎 耳 草 目 系 统 发 育 树 ， 涵 盖 虎 耳 草 目 下 全 部 
15 个 科 ， 各 科 均 具有 较 高 覆盖 率 ( 表 1) 。 
1.2 科 的 物种 多 样 性 
科 的 物种 多 样 性 数目 按照 生物 物种 名 录 数 据 库 (COL，https:Wwww.catalogueoflife.org/， 
访问 时 间 : 2021 年 5 月 ) 中 每 个 科 接 受 的 物种 名 称 进行 统计 。 变 种 和 亚 种 等 种 下 等 级 不 统计 。 
杂种 不 统计 。 由 于 物种 多 样 性 数目 不 符合 正 态 分 布 ， 参 考 Rabosky 等 (2012) 的 分 析 方 法 对 
科 的 物种 多 样 性 进行 了 以 10 为 底 的 对 数 转换 以 改善 拟 合 效果 。 
1.3 科 的 干 龄 和 冠 龄 计算 
冠 龄 (crown age) 与 干 龄 (stem age) 是 针对 系统 发 育 树 中 分 支 年 龄 的 两 个 不 同 的 概念 (图 1)。 
前 者 代表 该 类 群 现存 物种 最 近 共 同 祖先 的 年 龄 ， 而 后 者 代表 该 类 群 最 近 共 同 祖 先 与 其 姊妹 
类 和 群 的 最 近 共 同 祖先 的 年 龄 。 有 学 者 研究 认为 在 取样 较 好 的 情况 下 ， 冠 龄 用 来 做 数据 统计 
分 析 更 为 合理 ， 因 为 冠 龄 可 以 消除 由 于 进化 停滞 或 者 灭绝 事件 导致 的 长 枝 的 影响 (Sanchez- 
© Reyes et al., 2017)。 为 了 更 加 全 面 充分 地 评估 进化 时 间 与 物种 丰富 度 的 关系 ， 本 研究 同时 分 
析 了 各 个 科 干 龄 以 及 冠 龄 与 其 物种 多 样 性 的 关系 。 
为 了 提取 不 同 科 的 冠 龄 和 干 龄 ， 本 研究 利用 R 语言 ape 程 辑 包 (R Core Team, 2021; 
Paradis & Schliep, 2019) 中 getMRCAO 函 数 从 系统 发 育 树 上 找到 每 个 科 内 的 物种 对 应 的 共同 
祖先 的 节点 ， 该 节点 的 年 龄 即 为 该 科 的 冠 龄 ， 该 节点 的 父 节 点 即 为 该 科 的 干 龄 。 虎 耳 草 目 
个 节点 的 年 龄 使 用 R 程序 中 的 ape 程 辑 包 中 的 branching.times ( ) 函 数 提取 。 各 个 科 的 冠 
龄 节点 的 编号 使 用 ape 程 辑 包 中 的 getMRCAO 函 数 提取 。 根 据 每 个 科 的 冠 龄 节点 和 干 龄 节 
点 编号 (nodelable)， 分 别 获取 到 15 个 科 的 冠 龄 和 干 龄 。 
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冠 龄 Crown age (Myr) 


干 龄 Stem age (Myr) 


图 1 干 龄 和 冠 龄 的 示意 图 (单位 : 百 万 年 ) 
Fig. 1 Schematic diagram of “stem age” and “crown age” (Unit: Myr) 
1.4 科 的 多 样 化 速率 
本 研究 利用 宏观 进化 贝 叶 斯 分 析 法 (Bayesian Analysis of Macroevolutionary Mixtures, 
BAMMD) 使 用 及 语言 中 的 程 辑 包 BAMMtools 获取 (Rabosky et al., 2014)。 马 尔 可 夫 链 蒙特 卡 


> | 

») 
IA 
A J 
end 
T 


© 洛 (Markov chain Monte Carlo，MCMO) 运 行 长 度 为 1 000 万代， 每 1 000 代 取 样 一 次 。 运 行 
LO) 结束 后 ， 在 R 语言 中 使 用 程 辑 包 coda 检查 链 的 收敛 性 和 有 效 样本 (effective sampling sizes, 
e 

N ESS) 的 大 小 。ESS 检查 结果 大 于 200, WH MCMC 过 程 有 足够 的 独立 样本 ， 估 计 的 结果 比 


较 稳 定 。 

虎 耳 草 目 15 个 科 的 物种 多 样 化 速率 使 用 BAMMtools 程 辑 包 的 getCladeRates( ) 函 数 提 
取 。 该 方法 提取 了 科 内 所 有 分 支 的 枝 长 加 权 平 均 的 物种 形成 速率 (speciatiom) 和 灭绝 速率 

(extinction)。 净 多 样 化 速率 等 于 种 化 速率 与 灭绝 速率 的 差 值 。 由 于 灭绝 速率 的 估计 存在 较 

= 大 的 不 确定 性 ， 本 研究 中 同时 分 析 了 种 化 速率 和 净 多 样 化 速率 与 物种 多 样 性 的 关系 。 
© 此 外 ， 根 据 各 个 科 是 否 为 温带 适应 类 群 ， 利 用 ape 程 辑 包 中 的 drop.tip0 函 数 将 虎 耳 草 
目 系 统 发 育 树 拆 分 为 两 部 分 ， 即 温带 适应 类 群 和 常 绿 乔 木 类 群 ( 常 绿 乔 木 类 群 为 围 盘 树 科 
(Peridiscaceae)、 虎 皮 楠 科 (Daphniphyllaceae)、 鼠 刺 科 (Iteaceae) 和 曹 树 科 (Altingiaceae); 虎 
耳 草 目 下 其 余 11 个 科 为 温带 适应 类 群 )， 利 用 BAMMtools 程 辑 包 中 的 plot.bammdata() FK 2 
绘制 出 净 多 样 化 速率 随时 间 的 变化 图 像 。 
1.5 物种 多 样 性 与 多 样 化 速率 和 进化 时 间 的 关系 

本 研究 使 用 a Ame mnie ee 
型 (PGLS) 对 物种 丰富 度 与 多 样 化 速率 和 进化 时 间 进 行 了 分 析 。 为 验证 多 样 化 速率 假说 ， 
文 同 时 建立 了 种 化 速率 和 净 多 样 化 速率 与 物种 多 样 性 的 一 元 线性 回归 模型 和 PGLS 模型 。 
同时 也 分 别 建立 了 冠 龄 和 干 龄 与 物种 多 样 性 的 一 元 线性 回归 关系 和 PGLS 模型 来 验证 进化 
时 间 假 说 。 由 于 多 样 化 速率 和 进化 时 间 可 能 同时 对 物种 多 样 性 造成 影响 ， 因 此 ， 本 研究 也 
考虑 了 多 样 化 速率 和 进化 时 间 的 交互 效应 对 物种 多 样 性 的 影响 。 交 互 效 应 中 ， 进 化 时 间 由 
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冠 龄 来 代表 。LM 使 
gls0 函 数 完成 。PGLS 分 析 中 使 


个 科 只 有 1 个 代表 物种 。 
2 结果 与 分 析 
2.1 虎 耳 草 目 科 的 干 龄 和 冠 龄 与 物种 丰富 度 的 关系 
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FCMyr)， 而 最 古老 的 围 盘 树 科 (Peridiscaceae) 的 冠 龄 则 超 
EG 1)。 这 充分 体现 了 虎 耳 草 目 内 部 不 同 科 在 进化 时 间 
表 1 虎 耳 草 目 


] R 语言 中 的 其 而 


1 了 一 个 虎 耳 草 目 


15 个 科 有 关 数 所 


Sit 


函数 Im ( ) 完 成 。PGLS 使 用 R 程 辑 包 中 的 nlme 的 
的 科 级 系统 发 育 树 。 该 系统 发 育 树 是 使 用 


R 语言 的 “ape” 程 辑 包 中 的 drop.tip0 函 数 将 物种 数目 超过 1 个 科 中 的 物种 进行 去 除 ， 保 证 每 


E 轻 的 连 香 树 科 


尺度 上 的 差异 性 ( 表 D. 


Table 1 Relevant data statistics of 15 families of the order Saxifragales 


iiit TÉ 
Crown Stem age 
age (Myr) | (Myr) 
13.906 27.811 
41.726 83.452 
9.990 83.147 
3.009 66.357 
8.513 72.090 
82.662 110.769 
22.423 56.969 
63.851 91.030 
19.918 95.930 
12.216 108.265 
38.691 56.969 
28.338 72.090 
20.691 85.690 


总 物种 


Species 


number 


1 


104 


165 


195 


种 化 速率 灭绝 速率 
Speciation Extinction 
Tate (Myr!) rate (Myr!) 
0.053 0.011 
0.179 0.061 
0.052 0.020 
0.092 0.042 
0.051 0.020 
0.049 0.011 
0.172 0.092 
0.066 0.005 
0.079 0.021 
0.248 0.059 
0.218 0.119 
0.232 0.078 
0.232 0.041 


净 多 样 化 速率 


Net diversification 


Tate (Myr!) 
0.042 
0.118 
0.032 
0.050 
0.031 
0.038 
0.080 
0.061 
0.058 
0.189 
0.099 
0.154 


0.191 


取样 率 
Sampling 


fraction 
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REEF 64.599 85.690 869 0.128 0.013 0.115 0.66 
Saxifragaceae 
县 . 
n 79.134 106.364 — 1678 0.158 0.036 0.122 0.38 
rassulaceae 

M 分 析 的 结果 表明 冠 龄 对 物种 丰富 度 的 解释 率 可 以 达到 26.1%， 但 是 干 龄 对 物种 丰富 


度 的 解释 率 只 有 16.4%( 


0.01$5， 略 低 于 干 龄 与 物种 丰富 度 的 斜率 (05=0.019)。 但 无 论 是 冠 


影响 都 不 显著 (P>0.05)。 该 结果 说 明 虎 耳 草 目 内 15 个 科 的 物种 丰富 度 与 科 的 冠 龄 和 干 龄 没 


有 显著 相关 性 ， 进 化 时 


图 2). PGLS 分 析 的 结果 与 LM 一 致 ， 


冠 龄 与 物种 丰富 度 的 斜率 为 
了 龄 还 是 干 龄 对 物种 丰富 度 的 


间 不 是 影响 科 之 间 的 物种 多 样 性 差异 的 主要 因素 。 
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图 中 右 下 角 为 LM 模型 的 解释 率 (R”) 和 显著 性 检验 的 P 值 。 
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In the lower right corner are the LM model’s explanatory rate (R°) and P-value of significance test. 


图 2 虎 儿 草 目 15 个 科 的 冠 


二 乘 模型 (PGLS) 与 线性 回归 模型 (LMD) 拟 合 


Fig. 2 phylogenetic generalized least squares and linear regression model for crown age and 


R= 0.164 
P=0.134 


100 120 


里 龄 与 物种 丰富 度 (a) 及 干 龄 与 物种 丰富 度 (b) 的 系统 发 育 广 义 最 小 


species richness (a) and stem age and species richness (b) of 15 families within Saxifragales 


2.2 虎 耳 草 目 科 的 物种 多 样 化 速率 与 物种 丰富 度 的 关系 


虎 耳 草 目 净 多 样 化 速率 较 高 的 科 以 温带 和 高 山 类 群 为 主 ， 
， 例 如 茶 蕊 子 科 (Grossulariaceae)、 小 二 仙 草 科 (Haloragaceae) 和 景 天 科 
(Crassulaceae) 等 科 的 物种 形成 速率 和 净 多 样 化 速率 都 比较 高 ( 表 1, 
年 以 来 增长 较为 明显 (图 3， 图 人)。 物 种 主要 为 常 绿 乔 木 的 科 的 物 和 


物 为 主 的 草本 分 支 
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图 3)， 并 且 以 1 500 万 
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绿色 标注 的 科 是 常 绿 乔木 类 群 。 
The families shown in green are evergreen trees. 


3 虎 耳 草 目 净 多 样 化 速率 -时 间 图 


Fig. 3 Net diversification rate-time map of Saxifragales 
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Fig. 4 Net diversification rate curves over time for evergreen trees and temperate adaptive groups 


in the order Saxifragales 


LM 和 PGLS 分 析 发 现 ， 虎 耳 草 目 科 的 物种 形成 速率 和 净 多 样 化 速率 都 与 物种 丰富 度 


正 相 关 (图 5)。 但 物种 形成 速率 与 物种 多 样 性 


的 关系 在 LM 中 显著 (R=0.280，p=7.319,，P< 


0.05)， 而 在 PGLS 模型 中 接近 显著 (p=6.425，P<0.05)， 


0.05)。 


而 这 两 种 方法 都 发 现 净 多 样 化 速率 
与 物种 丰富 度 均 有 显著 正 相 关 性 (PGLS: 8=9.790，P=0.031; 


LM: R=0.380, p=11.791, P< 
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阴影 代表 LM 模型 的 95% 的 置信 区 间 。 图 中 右 下 角 为 LM 模型 的 解释 率 (R") 和 显著 性 检验 的 P 值 。 下 同 。 


The shadow represents 95% confidence interval of the LM-model. In the lower right corner are the LM model’s 


explanatory rate (R?) and P value of significance test. The same below. 
图 5 虎 儿 草 目 15 个 科 的 种 化 速率 (a) 和 净 多 样 化 速率 (b) 与 物种 丰富 度 的 系统 发 育 广义 最 小 
二 乘 模型 (PGLS) 与 线性 回归 模型 LMD) 拟 合 
Fig. 5 Phylogenetic generalized least squares model (PGLS) and linear regression model (LM) 


fitted speciation rates (a) and net diversification rate (b) with the species richness of 15 families in 


the Order Saxifragales 


2.3 进化 时 间 与 多 样 化 速率 的 综合 效应 与 科 物 种 丰富 度 的 关系 
进化 时 间 与 净 多 样 化 速率 对 物种 丰富 度 的 共同 作用 在 LM 和 PGLS 模型 中 都 显著 (LM 


R’ =0.603, f=0.297, P«0.05; PGLS: 6=0.287，P < 0.05) (图 6). LM 模型 中 ， 进 化 时 间 和 
进化 速率 的 共同 作用 解释 了 物种 多 样 性 变化 的 60%，PGLS 模型 的 斜率 也 比 单独 使 用 进化 
时 间 和 多 样 化 速率 高 。 该 结果 支持 了 进化 时 间 与 多 样 化 速率 的 综合 效应 对 于 虎 耳 草 目 科 间 
物种 丰富 度 差异 的 影响 ， 且 比 多 样 化 速率 这 一 单一 变量 对 其 的 影响 的 解释 率 更 高 。 


线性 回归 模型 LM A 


e 
— 系统 发 育 广 义 最 小 二 乘 模型 PGLS Ex 
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物种 丰富 度 (对 数 ) Species richness (Log) 
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冠 龄 与 净 多 样 化 速率 的 乘积 


Crown age (Myr) * Net diversification rate (Myr') 
图 6 虎 儿 草 目 15 个 科 进 化 时 间 和 净 多 样 化 速率 的 综合 效应 与 物种 丰富 度 的 系统 发 育 广 义 最 
小 二 科 模 型 (PGLS) 与 线性 回归 模型 (LM) 拟 合 
Fig. 6 Phylogenetic generalized least squares model (PGLS) and linear regression model (LM) 
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fitted the interaction effects of evolutionary time and net diversification rate with the species 


richness of 15 families in the Order Saxifragales 


3 讨论 与 结论 T5 
在 本 研究 中 ，PGLS 分 析 发 现 虎 耳 草 目 下 15 个 科 的 物种 多 样 性 与 净 多 样 化 速率 呈现 显 


著 的 正 相 关 关 系 ， 但 是 与 科 的 冠 龄 正 相 关 关 系 不 显著 。 线 性 回归 模型 中 ， 净 多 样 化 速率 对 
物种 多 样 性 的 解释 率 达 到 38%。 由 于 Scholl & Wiens (2016) 发 现 多 样 化 速率 与 年 龄 的 交互 作 
用 可 能 是 导致 物种 多 样 性 差异 原因 ， 因 此 本 研究 也 检验 了 物种 净 多 样 化 速率 与 冠 龄 的 交互 
作用 ， 并 发 现 交互 作用 对 物种 多 样 性 的 解释 率 达 到 60%。 因 此 本 研究 认为 净 多 样 化 速率 的 
差异 是 主导 物种 多 样 性 差异 的 主要 因素 ， 进 化 时 间 与 多 样 化 速率 的 共同 作用 在 虎 耳 草 目 物 
种 多 样 性 差异 中 起 到 了 次 要 作用 。 

本 研究 的 结果 不 支持 进化 时 间 假 说 ， 这 与 Rabosky 等 (2012) 和 Scholl & Wiens (2016) 得 
到 的 结论 一 致 ， 即 不 同类 群 间 物 种 丰富 度 的 差异 并 非 是 由 进化 时 间 的 差异 引起 。 尽 管 
Pyron & Wiens (2013) 的 研究 中 发 现 热 带 两 栖 类 物种 类 群 之 间 的 物种 多 样 性 差异 与 冠 龄 有 显 
著 的 正 相 关 关系 ， 但 是 在 该 研究 中 也 同时 发 现 物种 多 样 性 与 净 多 样 化 速率 有 着 显著 的 正 相 
关 关 系 。 该 研究 还 发 现 物种 多 样 化 速率 与 净 第 一 性 生产 力 和 面积 存在 正 相 关 关 系 。 热 带 地 
区 的 科 的 物种 的 生态 位 演化 快 ， 种 化 速率 高 ， 灭 绝 速率 低 ， 古 老 类 群 能 够 保存 ， 是 该 热带 
类 群 物种 多 样 性 高 的 重要 原因 。 在 虎 耳 草 目 中 ， 草 本 和 灌木 等 能 够 适应 温带 环境 的 科 有 着 
更 高 的 物种 形成 速率 和 净 多 样 化 速率 ， 数 据 分 析 也 显示 这 可 能 是 虎 耳 草 目 中 物种 丰富 度 与 
J 净 多 样 化 速率 正 相 关 关 系 形成 的 主要 原因 。Folk 等 (2019) 对 虎 耳 草 目 生态 位 祖先 重建 的 结 
> 果 表明 景 天 科 、 虎 耳 草 科 和 茶 万 子 科 等 温带 适应 的 科 的 共同 祖先 ， 在 8 千 万 年 前 就 已 经 能 
c 够 适应 较为 干旱 和 寒冷 的 环境 。 在 渐 新 世 全 球 气候 变 冷 时 ， 这 些 类 群 的 物种 多 样 化 速率 就 
已 经 高 于 虎 耳 草 目 木 本 分 支 的 曹 树 科 等 。 中 新 世 后 ， 随 着 全 球 气 候 的 进一步 变 冷 ， 这 些 温 
带 类 群 的 物种 多 样 化 速率 快速 上 升 。 相 比 喜 湿热 环境 的 草 树 科 、 虎 皮 楠 科 等 类 群 ， 虎 耳 草 
» 目 中 的 草本 和 灌木 类 群 对 寒冷 环境 的 良好 适应 可 能 为 该 类 群 在 全 球 变 冷 背景 下 的 快速 多 样 
化 提供 了 契机 。 例 如 ， 小 二 仙 草 科 (Haloragaceae) 在 始 新 世 - 渐 新 世 的 全 球 气候 变 冷 和 和 干旱 化 

背景 下 ， 分 布 范 围 不 断 扩张 ， 多 样 化 速率 经 历 了 大 幅度 提升 (Moody & Garcia, 2021). 
地 质 历 史 事件 不 仅 会 重新 塑造 原 有 的 地 形 地 貌 还 会 改变 区 域 气候 条 件 ， 从 而 对 其 间 的 
生物 演化 过 程 产生 影响 ( 唐 志 尧 和 方 精 云 ，2004; Antonelli et al, 2018)。 山 脉 的 隆 升 直 接 导 
致 的 种 群 隔离 分 化 以 及 间接 通过 生境 的 复杂 化 等 导致 的 种 群 的 生殖 隔离 和 丰富 的 生态 位 促 
进 了 高 山 植物 的 快速 演化 ， 这 也 与 虎 耳 草 目 多 个 科 的 快速 物种 多 样 化 有 着 密切 的 联系 
(Zhang et al., 2014a; 刘 杰 等 , 2017; Moody & Garcia, 2021)。 山 地 草本 类 群 对 寒冷 生境 的 快速 
适应 和 分 化 ， 也 是 青藏 高 原 、 北 美 西海 岸 山脉 和 欧洲 南部 阿尔 卑 斯 山脉 成 为 虎 耳 草 目 中 高 
山 类 群 物种 多 样 性 中 心 的 重要 原因 之 一 (Zhang et al., 2008; Zhang et al., 2014b; Ebersbach et 
al., 2017; Xing et al., 2017; Stubbs et al., 2020)。 随 着 青藏 高 原 的 隆 升 ， 许 多 植物 类 群 出 现 了 
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快速 的 多 样 化， 其 中 包括 虎 耳 草 目 中 适应 寒冷 的 类 群 ， 虎 耳 草 属 (Saxifraga L.)、 红 景 天 属 
(Rhodiola L.) 等 类 群 的 多 样 化 速率 快速 升 高 (Zhang et al., 2014a; Ebersbach et al., 2017; Folk et 
al, 2019)。 因 此 ， 虎 耳 草 目的 高 山 类 群 物种 多 样 性 快速 增加 ， 物 种 多 样 性 与 净 多 样 化 速率 
呈现 正 相 关 的 关系 。 

对 类 蓄 项 目 Rosids) 的 研究 发 现 ， 包 括 虎 耳 草 目 在 内 的 其 他 目的 多 样 化 速率 也 随 着 全 球 

变 冷 而 上 升 ， 并 且 多 样 化 速率 随 着 纬度 升 高 而 升 高 (Sun et al., 2020). Igea 等 (2020) 针 对 被 子 
植物 的 研究 也 发 现 类 似 的 现象 ， 即 植物 类 群 在 温带 和 寒带 的 多 样 化 速率 常常 高 于 热带 。 这 
些 研 究 的 结果 表明 温带 类 群 在 全 球 变 冷 的 过 程 中 物种 多 样 性 在 不 断 地 积累 。 尽 管 空间 上 看 ， 
由 于 高 纬度 受 冰 期 影响 严重 ， 物 种 多 样 性 较 低 ， 但 是 从 类 群 上 看 ， 温 带 类 群 的 物种 多 样 性 
的 增加 速率 高 于 热带 类 群 。 因 此 ， 除 了 虎 耳 草 目 外 ， 在 被 子 植物 的 科 中 也 可 能 存在 类 似 的 
现象 ， 即 科 间 的 物种 多 样 性 差异 与 净 多 样 化 速率 可 能 存在 显著 的 正 相 关 关 系 。 
有 研究 者 提出 物种 多 样 性 差异 产生 的 真实 原因 很 可 能 是 多 因素 共同 导致 的 ， 分 支 年 龄 
和 多 样 化 速率 可 能 都 在 其 中 发 挥 着 作用 (Valente et al., 2011; Bloom et al., 2014; Yan et al., 
2018)。 近 年 来 ， 也 有 学 者 提出 该 问题 的 讨论 应 当 换 一 种 思考 模式 ， 由 验证 哪 种 假说 更 正确 
转向 研究 这 些 影 响 因 素 之 间 的 互 作 关系 ， 并 建立 相关 模型 来 更 加 数据 化 与 统计 学 化 地 看 待 
物种 丰富 度 差异 的 形成 机 制 (Pontarp et al., 2019)。 尽 管 本 研究 中 发 现 了 净 多 样 化 速率 与 年 龄 
的 交互 作用 对 物种 多 样 性 差异 有 很 高 的 解释 率 ， 但 是 目前 该 结论 是 否 适用 于 其 他 类 群 ， 还 
有 待 验 证 。 总 之 ， 迄 今 为 止 该 问题 仍然 留 有 许多 值得 讨论 的 空间 ， 利 用 更 为 复杂 的 统计 模 
型 方法 以 及 更 为 全 面 的 动 植物 数据 集 去 探讨 该 问题 可 能 是 未 来 科研 工作 者 们 的 研究 方向 。 

核 基因 的 同 义 突变 率 比 叶绿体 基因 高 5 倍 ， 比 线粒体 基因 高 20 倍 ， 且 核 基 因 拥 有 多 个 
独立 位 点 、 双 末 遗 传 等 特征 (Small et al., 2004)， 因 此 ， 采 用 核 基 因 构 建 的 系统 发 育 树 与 以 
往 小 片段 建树 相 比 ， 能 够 更 加 全 面 和 客观 地 反映 分 支 的 系统 发 育 关系 。 本 研究 中 使 用 的 系 
= 统 发 育 树 是 Folk 等 (2019) 基 于 301 MAF MEM EA BB 3E EA) SEE ASK BR AT SS 

时 间 的 基础 上 ， 结 合 叶 绿 体 基因 组 小 片段 和 ITS 序列 所 构建 。 该 研究 结果 中 各 个 科 的 关系 
与 APG IV (Chase et al., 2016) 一 致 。 然 而 ， 由 于 物种 采样 困难 以 及 大 量 的 核 基 因 测序 时 间 
和 经 费 也 比较 多 ， 因 此 核 基因 组 或 者 大 量 核 基因 建树 的 物种 覆盖 率 都 比较 低 。 而 取样 率 低 
对 物种 多 样 化 过 程 的 分 析 的 准确 性 有 着 非常 大 的 影响 (Chang et al., 2020)， 因 此 基因 组 结合 
小 片段 建树 是 提高 物种 覆盖 率 的 解决 方案 之 一 。 然 而 ， 由 于 小 片段 的 信息 位 点 不 足 ， 可 能 
不 能 够 十 分 准确 地 估计 近 缘 物种 之 间 的 系统 发 育 关 系 ， 也 可 能 对 系统 发 育 树 的 枝 长 估计 产 

影响 ， 从 而 影响 到 多 样 化 速率 的 估计 。 随 着 叶绿体 基因 组 测序 价格 的 下 降 和 分 析 方 法 的 
成 熟 ， 增 加 叶绿体 基因 组 的 物种 覆盖 率 可 能 是 解决 近 缘 物种 系统 发 育 关系 的 较为 经 济 快捷 
的 方法 。 
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